Optimal Dynamical Decoupling Sequence for Ohmic Spectrum 
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We investigate the optimal dynamical decoupling sequence for a qubit coupled to an ohmic envi- 
ronment. By analytically computing the derivatives of the decoherence function, the optimal pulse 
locations are found to satisfy a set of non-linear equations which can be easily solved. These equa- 
tions incorporates the environment information such as high-energy (UV) cutoff frequency w C) giving 
a complete description of the decoupling process. The solutions explain previous experimental and 
theoretical results of locally optimized dynamical decoupling (LODD) sequence in high-frequency- 
dominated environment, which were obtained by purely numerical computation and experimental 
feedback. As shown in numerical comparison, these solutions outperform the Uhrig dynamical 
decoupling (UDD) sequence by one or more orders of magnitude in the ohmic case. 
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I. INTRODUCTION 



Suppressing decoherence is one of the fundamental is- 
sues in the field of quantum information processing. De- 
coherence, which has been caused by the environmental 
noise, plagues almost all the implementations of quantum 
bit. To eliminate the unwanted coupling between a qubit 
and its environment, several schemes have been proposed 
and tested. Among them a promising one is dynamical 
decoupling [lj-0] , which restores the qubit coherence by 
applying delicately designed sequence of control pulses. 

For a qubit that can be modeled by a spin-1/2 par- 
ticle, the oldest dynamical decoupling sequence is peri- 
odic dynamical decoupling (PDD). Originated from pulse 
sequences widely used in nuclear magnetic resonance 
(NMR) [|, the PDD sequence consists of periodic and 
equidistant ir pulses. To achieve better performance, 
there has been an extensive study in how to optimize 
the pulse locations [r34IH . One important progress is the 
powerful Uhrig DD (UDD) [l(J, which employs n pulses 
located at tj according to the simple rules 

^• = sin 2 (j7r/2(n + l)), 

where 6j — tj/T and T is the total evolution time. UDD 
is first derived on spin-boson model and further proved 
to be universal in the sense that it can remove the qubit- 
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bath coupling to n order in generic environment 

Beyond UDD, another locally optimized dynamical de- 
coupling (LODD) sequence has drawn great attention 
[l4j. LODD, along with its simplified version optimized 
noise-filtration dynamic decoupling (OFDD) [lfl, gener- 
ates the decoupling sequence by directly optimizing the 
decoherence function using numeric methods as well as 
experimental feedback. It has been shown to be able to 
suppress decoherence effect by orders of magnitude over 



UDD for certain noise spectrum, especially for the one 
with a high frequency part and sharp high-energy (UV) 
cutoff. 

However, in spite of the great experimental success, an- 
alytical results about the LODD sequence is insufficient. 
Until recently S. Pasini and G. S. Uhrig has made an 
analytical progress in optimizing the decoherence func- 
tion for power law spectrum (PLODD) The power 
law spectrum uj a for a < 1 without UV cutoff is consid- 
ered. They minimize the decoherence function through 
expanding the function and separating, canceling diver- 
gences from the relevant terms and solving variation 
problems. Inspired by Pasini's work, we try to analyze 
the LODD problem with respect to the ohmic spectrum 
S(lu) ~ to and a sharp UV cutoff. Ohmic noise is the ma- 
jor decoherence source often found in a qubit's environ- 
ment, for example, the semiconducting quantum dot [l6j 
and superconducting qubit j 1 71 ] - Optimal performance 
pulse sequence is found analytically which entirely dif- 
fers from the UDD sequence in such environment. Wc 
call this kind of optimal sequence HLODD (LODD for 
ohmic spectrum) for short. 

We organize this paper as follows. In the second sec- 
tion we propose the optimization problem of the deco- 
herence function. In section III, we derive the analytical 
equations for the optimal pulse sequence. In the follow- 
ing section, we run a simulation to verify our results. 
Conclusions are put in section V. 



II. OPTIMIZATION OF THE DECOHERENCE 
FUNCTION 



Given a two-level quantum system, when the environ- 
mental noise behaves quantum-mechanically, we use the 
long-established spin-boson model with pure dephasing 
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Here we ignore the qubit free evolution hamiltonian. On 
the other hand, when the qubit is subjected to classical 
noise, the system is modeled as [H, fl9| 



H = -flu!- 
2 l 



■P(t)]<r» 



(2) 



where the O is the qubit energy splitting and /3(f) the 
classical random noise. Let t be the total evolution time, 
and n pulses are applied at ti < t 2 < ... < t n in sequence 
with negligible pulse durations . We use the notation 
Sj = — . This naturally leads to the definition of to = 
and t n +i = 1. In either (1) or (2), the decay of coher- 
ence under the dynamical decoupling sequence can be 

described by the decoherence function [Io|, 0, Ho|, HJ 
e -2 X (t) witn 



X(t)= I™ ^\y n (ut)\ 2 du, 
Jo w 



(3) 



where S(uj) is environmental noise spectrum. The filter 
function y n (t) is given by 



y n (t) = ! + (-!) 



2^(-l)V"^. (4) 



3=1 



Thus minimization of with respect to Sj gives the 
optimal decoupling sequence. 

We now consider the case when the noise spectrum is 
ohmic with a sharp cutoff at lj c , i.e. S(ui) = Sou;0(w c — 
to). So is an irrelevant constant factor and is unit step 
function. Then minimization of (3) turns to minimization 
of /„ with 



\Vn{z)[ 



dz, 



(5) 



where z c = u) c t. Since y n (0) — 0, the IR convergence 
insures the integral converges to a finite value [HI]. 



III. DERIVATION OF OPTIMAL PULSE 
SEQUENCE 

We follow the approach of Pasini and Uhrig [l3j to 
treat the integral (5). Here we use notation 

_ / if j = 0,n + l, 
qj ~\ 1 ifje{l,2,...,n}, 



and 



from which we get 



Ay = i(Si -Sj), 



n+l 



\y n {z)\ 2 = £ 2*+*(-l)*+'e* A «. 

i,j=0 



Then the integral /„ can be expressed as 
I n = lim I n (x), 



n+l 



I n (x)= J2 2^(-l) i+ %(x), 

i,j=0 



(6) 



where the integrals Iij{x) are 



Iij(x) = 



-dz 



-A tJ z c e - z 



AijX 



(7) 



The limit x — > + is carried out because Iij (0) does not 
exist for arbitrary Making use of the series represen- 
tation of exponential function [22j 



Ei{z) 



-dt 



= -7 ' 



k=l 



(~l) fc+1 k 

m z ' 



where 7 is the Euler-Mascheroni constant and the sum 
converges for all the complex z, hj{x) can be written as 

Iij(x) = Ei(-Aijx) - Ei(-Aijz c ) 
°° A fc 

= l n (z c /x) + -x k ). (8) 

fe=i 

Since we always have j/«(0) = which implies 

n+l 

\y n {0)\ 2 = ]T 2*+*(-l) i+ ' =0, 
we can now proceed by taking the limit ir — > 0+ in 7„ 



(9) 



7 « = E 2^(-ir + - ? '[ln(z c / 2 ;)+^ I |(z c fe - a; fc )] 

^ i,j=0 fe=l 
n+l 00 A fe 

= iim y 2*+* y ^(j** - 



«,i=o 

71 + 1 OO 

ij=0fc=l ' 



fc=l 

A? 



(10) 



To minimize I n , UDD requires the first n derivatives of y n 
vanish while OFDD simplifies the optimization process 
by replacing S(uj) by a constant. Here we attempt to 
minimize I n directly to obtain optimal pulse sequence. 
We notice that at the optimal pulse locations Sj (j — 
l,2,...,n), the gradient of /„ vanishes. So we impose 



the following conditions 



ai„ 



0, for m from 1 to n. 



Although (10) are complex infinite series, we can still 
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explicitly compute the derivatives of (10) as long as these 
derivatives converge. For arbitrary m we have 



FIT 8 n+ °° A k - 



i,j=0 k=l 
n+1 oo 



klk 



i 

n+1 

+EE* 



A fc 

q m +qi ( i \m+i ^mt fc\ 

1 j fc!fc Zcl 



i=0 fc=l 
n+1 oo 



i=0 fe=l 



(ii) 



The terms with k odd cancel, so the result can be sim- 
plified as 



n+1 oo 2k 
^^ 29m + 9 , + l ( _ ir+i ^. i 2 fc((5m _^ ) 



i=0 k=l 
n+1 



(2fc) 



\2fe-l 



E 

n+1 



-2«»+«+ l (-i) wl+i E 7#v i2fc(,5m - ^ 



x2/c 



V F ^2'"»+«+ 1 (-l) m+i {cos[(5 m - fc)z c ] - 1}. 



(12) 



Here we have used the expansion cos (z) = 



k=0 



(-1)* -2fc 
(2*:)! Z 



which converges on the whole complex plane. From (12) 
we know that the derivatives of I n indeed converge to a fi- 
nite value. Thus the optimal pulse locations {5\, 62, ■■■8 n \ 
shall satisfy the following non-linear equations 



n+1 



J2 ir —ir2 qm+q > +1 (-iy n+l {cos[(S m - 6 t )z c ] - 1} = 0. 



(13) 

Equations (13) are main results of this paper. The opti- 
mal sequence obtained from (13) is quite different from 
the UDD sequence obeying 




0.9 - 
0.8 - 
0.7 - 
0.6 - 
0.5 - 
0.4 - 
0.3 - 
0.2 - 
0.1 - 



1 1— 
0.9 - 
0.8 - 
0.7 - 
0.6 - 
0.5 - 
0.4 - 
0.3 - 
0.2 - 
0.1 - 









FIG. 1: Comparison between UDD and HLODD sequence for 
different UV cutoff frequency oj c . Pulse sequences Si for n — 2 
and n — 5 are plotted in one figure under the same w c . 



n+1 

^2*(-l)^ = 

J=l 

for p — {1, 2, ...n}. For the ohmic spectrum, our equa- 
tions incorporate the UV cutoff frequency uj c , indicating 
that the solutions are specially tailored to combat this 
kind of noise. Although the UDD sequence is universal 
in suppressing decoherence, we believe that the HLODD 
sequence will outperform the UDD sequence in the ohmic 
environment. In the next section, we use numeric meth- 
ods to illustrate the performance of HLODD sequence. 



IV. NUMERICAL RESULTS 

We start our simulation by solving the non-linear equa- 
tions (13) and evaluate the decoherence function with 
these solutions. First, we set the total evolution time 
t = 1 without loss of generality. Then z c = lu c and we 
can concentrate on analyzing the influence of the cutoff 
frequency ui c . Computing solutions to (13) for differ- 
ent w c , we find that the optimal pulse sequences behave 
differently. We also evaluate the UDD sequence for com- 
parison. 
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FIG. 2: Semi-log plot of /„ versus pulse number n. UDD and 
HLODD sequences are compared. 

As shown in Fig. [T] deviation of the pulse locations 5, 
in HLODD sequence from their UDD counterparts in- 
creases with uj c . This agrees with our intuition since 
UDD focuses on suppressing decoherence by minimiz- 
ing |yn(z)| in the neighborhood of y n {0), weakening its 
ability to maintain small |?/„(^)| on the other end of the 
spectrum. For large oj c , UDD sequence is no longer op- 
timal. In addition, we can see pulse number n plays an 
important role. By increasing n, UDD can narrow the 
difference from HLODD. The difference between the two 
sequences when n = 2 is greatly reduced when we in- 
crease n to 5, see Fig. [TJ Especially for the case oj c = 1, 
the difference is completely removed. However, for larger 
to c this gap can't be removed by increasing n. 

Next, to demonstrate the optimal decoupling ability of 
HLODD sequence, we compute /„ versus n while uj c is 
chosen to be 5. The results are depicted in Fig. [2J and 
again are compared with UDD. The obtained solutions 
yield a significant improvement over UDD. For fixed n, 
the HLODD suppresses decoherence better than UDD 
by one or two orders of magnitude which is in agreement 
with the results in H^, where LODD and OFDD 
sequences are tested for 9 Be + qubits in a penning ion trap 
and various spectrum. The qubit error rates are below 
1CP 5 when n > 5, and we see that HLODD is capable of 
suppressing the error rates far below the Fault- Tolerance 
error threshold [23] by increasing n. Furthermore, by 



inspecting the points on the HLODD curve, we expect 
the HLODD sequence suppresses decoherence in power 
law n as UDD. 

At last, we would like to explain the numerical results 
in another way. If we fixed UV cutoff frequency uj c at 
the beginning, and compare the HLODD performance 
for t = 1, t = 5, and t = 10, the numerical results would 
be the same since z c did not change. So we can also 
conclude that for the same number of pulses n, HLODD 
will beat UDD with increasing total evolution time t. 
V. CONCLUSIONS 

In this paper we analytically find the optimal pulse 
locations to decouple a qubit in an ohmic environment. 
By deriving the analytical expressions for the derivatives 
of decoherence function, we obtain a set of non-linear 
equations which the optimal pulse sequence must obey. 
These equations are completely different from UDD and 
are more accurate, because they incorporate the effect of 
UV cutoff frequency uj c . 

In our numerical simulation, the analytical results pro- 
vide an improvement over UDD sequence by an order or 
two of magnitude, which is consistent with previous re- 
sults in LODD and OFDD obtained by purely numerical 
minimization and experimental feedback. We have to 
mention that the pulse performance is influenced by the 
sharp UV cutoff frequency w c greatly. The larger the UV 
cutoff av, the more HLODD deviates from UDD. Early 
work [l3|, [13, [H, H3] has pointed out that for soft large 
UV cutoff, UDD performs even worse and LODD is still 
a better choice. However, the integral (3) for S(u>) with 
a soft cutoff is hard to analyze. 

In conclusion, our work provides an analytical solution 
to optimal dynamical decoupling for ohmic case. Our 
derivation is based on ohmic spectrum, but we believe it 
can be extended to super-ohmic case S(u>) ~ ui a {a > 1) 
via slight modification. 
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